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> ' Abstract 

cn 

' We argue that Einstein gravity coupled to a Born-Infeld theory provides an attractive 

candidate to represent dark matter and dark energy. For cosmological models, the Born- 

■ Infeld field has an equation of state which interpolates between matter, to = (small 
00 ' 

«) I times), and a cosmological constant lo = —1 (large times). On galactic scales, the Born- 



Infeld theory predicts asymptotically flat rotation curves. 

1 Introduction and summary of results 

The existence of dark matter and dark energy is now firmly established phenomenologically 
P,[2] but the theoretical understanding is far from complete. Einstein equations require "exotic" 
components in the right hand side corresponding to about %96 of the total energy density today. 
Understanding the microscopic nature of these extra components is one of the most challenging 
and important problems faced by theoretical physics at present. 

Supersymmetric and exotic particles in the standard model are the best candidates for dark 
matter (for a review see [3j). Several experiments are now being devised for a direct detection 
of these particles. Alternative descriptions based on modifications to gravity have also been 
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explored with interesting results. See IH El El El El IS] , and references quoted therein, for some 
of these efforts and its consequences. For a recent review of the Einstein aether theory see |10j . 

The problem of dark energy is somehow more recent, although the issue of the cosmological 
constant has been around for a long time. The discovery of an accelerating Universe [11] 
resulted in deep changes in cosmology. The simplest explanation for this phenomena is a small 
positive cosmological constant, but many other possibihties have been explored (see [121 [13] for 
recent reviews). 

In this paper we consider an action for general relativity coupled to a Born-Infeld theory. 
The Born-Infeld theory has as fundamental variable a symmetric connection C^^^(x). C'^j^p has 
the same symmetries and transformation properties of the Christoffel symbol but is independent 
from it. The action is 

+ 1 d'xCm{'^,gp,), (1) 

where \Afj^^\, for any A^j^^, denotes the absolute value of the determinant of A^,^. K^^, is the 
"Ricci" curvature associated to C^^p(x), 

K,. ^ (K^.^p = + C^^aClp - [a - /3]). (2) 

Besides Newton's constant, the action ([1]) has two extra parameters: / is a length and a is 
dimensionless. \l/ denotes all baryonic fields and Cm the baryonic Lagrangian. 

The action ([T|) is similar in spirit although different in interpretation to the Born-Infeld 
gravity action proposed by Deser and Gibbons [H], 

I[9,u] = J ^J\9,u-l'R,u + X^,{R)\ (3) 

and elaborated in [15]. As discussed in [T3], the term X^y{R) must be chosen such that the 
action is free of ghost, and free of Schwarzschild-like singularities. The action ([3]) is an action for 
pure gravity, and it can be seen as a natural extension to spin two of the scalar \J\gfiu + d^(l)du4>\ 
and vector ^Jlg^u + F^^ul Born-Infeld (BI) theories. For the scalar and vector BI theories the 
equations of motion are of second order. For the spin two theory this is not automatic and 
requires the addition of X^^. 
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The action ([T]), on the other hand, gives rise to second order equations because K^^{C) 
depends on first derivatives of the field C^^^. This action, however, is not an action for pure 
gravity but gravity coupled to C^^^. The equations of motion are discussed in the appendix 
and in Sec. [2] below. 

It is known (e.g. [16]) that general relativity with cosmological constant is dual to Edding- 
ton's action [24] I[C] ~ / ^J\K^\. The action ([1]) can then be interpreted as general relativity 
interacting with its own dual field theory. 

The action ([1]) can also be motivated by looking at general relativity without metric |22j . 
This interpretation will be discussed in Sec. |5l 

Our main goal in this paper is to argue that the field C^^^ has good properties to represent 
dark matter and dark energy. We shall study the equations of motion following from ([1]) and 
prove the following properties. 

1. For a cosmological model, there exist solutions where the expansion factor a{t) behaves 
as a{t) ~ e^* for large t, and as a{t) ~ t^/^ for small t. The equation of state for the fluid 
interpolates between p = and p = —p. The parameters in the solution can be adjusted 
such that this field contributes to ~ 23% of the total matter energy density and ~ %73 
of vacuum energy density, as required by observations 0. 

2. For a spherically symmetric configurations, the action ([1]) predicts asymptotically fiat 
rotation curves, as required by galactic dynamics. The parameters involved in this solution 
can also be adjusted to deal with realistic situations. 

We would like to stress the simplicity of this proposal. The "Born-Infeld" term is all we need 
to account for both dark energy and dark matter, at least for the problems described above. 
More complicated tests, like lensing, fluctuations, and others will be discussed elsewhere [191 120] • 
See also [21]. 

^Couplings between dark matter and energy have appeared in [T7], and in [TB] involving a Chapligyn gas. 
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2 The Equations of Motion 



2.1 A bi- metric theory 

The fields varied in the action ([1]) are the metric g^i, and the connection C^^^p. Both fields 
are independent. At the level of the equations of motion, the connection C^^^ can be written 
in terms of a second metric q^^p. (The full action can also be written as a bi-metric theory 
[23].) This action then represent a bi-metric theory if gravity. This result follows closely the 
structure of Eddington's theory [2l]. We shall postpone a detailed derivation for the appendix 
and include here only the result. 

Let q^u{x) be a rank two invertible symmetric tensor satisfying the metricity condition 

Dpq^u = (4) 

with respect to C^'^p. Since C^^p is symmetric this implies C^^^ = \q^'^{qav,p + qap,y — Qu,p,a)i and 
for every q^^^ there is a unique C^^p. 

The equations of motion derived from the action ([1]) can be written completely in terms of 
Qpy and g^i,, and take the very simple form 

G,u = -L^g,aq'"'g(3u + 87iGr-;:, (5) 

Kpu = j^{9f,u + aqf,u) (6) 

T^pl is the energy momentum tensor associated to the baryonic Lagrangian q'^'^ is the 

inverse of The derivation of these equations is left for the appendix. 

Equation (jSj) is the Einstein equation. The first term in the right hand side is the contri- 
bution from the Born-Infeld action. Our main goal will be to prove that this fiuid can account 
for dark matter and dark energy. 

2.2 The de-Sitter solution 

The de-Sitter spacetime is an exact solution to this theory. This can be seen as follows. (The 
de-Sitter spacetime is expected to be relevant after matter becomes negligible so we set here 
T(-) = 0.) 



Suppose there exists solutions of the equations of motion with R^^, = Ag^^. It is direct to 
see that this imphes that both metrics must be proportional, 

Qi^uix) = g^„y{x) (7) 

with 7 a constant. The constant 7 can be computed as follows. Replacing in ([6]) we derive, 

Rfiu = ^ (7" + 1) 9f^u- (8) 
Replacing in (with T^;^) = 0) we derive 

Consistency determines 7, 

Thus, the Born-Infeld field can behave as a cosmological constant with the value 

The value a = 1 is a critical point where cosmological solutions ceases to exist. Curiously, we 
shall see that a good fit for the Friedman equation requires a to be close, but not equal, to one. 



3 Friedman cosmological models 

The evolution equation for the scale factor in fiat cosmological models is given by the Friedman 
equation (neglecting radiation) 



+ ^K. (12) 



Current values for the (relative) densities of barionic matter fi^m, dark matter Vtdm and vacuum 
energy are, 

VLb^ ^ 0.04, VLdm. ^ 0.23, ~ 0.73. (13) 



Among the components appearing in the right hand side of ( ]T2l) . only the ~ 0.04 fraction of 
baryonic matter is theoretically well-understood. The other 0.23 + 0.73 = 0.96 fraction remains 
a great mystery. 



3.1 Goal of this section 



The goal of this section is to demonstrate that the field C^^^ behaves like dark matter for small 
times, and as dark energy for larger times. In other words, its equation of state evolves from 
p = into p = —p. Adjusting the parameters a and /, plus initial conditions, the Born-Infeld 
field can account for both the Qdm and Qa contributions in (|T2l) . Thus, the action ([T]), is capable 
to reproduce the correct evolution of the scale factor without adding neither dark matter nor 
dark energy. 

Our approach does not shed any light into the particular values for ^a, ^dm, ^bm and other 
cosmological parameters. We shall only prove that / and a can be chosen such that the pre- 
dictions from ([1]) are consistent with the Friedman equation f|T2l) . In particular we have chosen 
here to set k = and consider only fiat models. There is no particular reason for the choice 
other than simplicity. A full analysis with a varying k and including other developments will 
be reported in |19j . 

3.2 The ansatz and equations 

To solve ([5]) and ([6]) we assume that both g^^, and are homogeneous, isotropic and with fiat 
spatial sections. Using the gauge freedom in the time coordinate to fix gu = — 1, the ansatz for 
g^„ and is then, 

g^.ydx^dx" = -dt^ + a{tf{dx^ + dy^ + dz"^), (14) 
q^^dx^'dx" = -X{tfde + Y(tf{dx^ + dy^ + dz^) (15) 

where a{t),X[t),Y{t) are arbitrary functions of time to be fixed by the equations of motion 
and initial conditions. 

As usual for fiat models, and to match the choice made in (fT2l). we set 



a(t)|todaH = 1, -f^O = O(i)|ioday (16) 

and use Hq to define a natural dimensionless time coordinate HqI. The time coordinate in all 
expressions from now on refer to this choice. 

Equations ( 15)161) for the ansatz fll4|15l) become, 

1 p 



a? 3PH^ X p. 



(17) 



?)XYaa 



plus second order equations related to (1171119^ by Bianchi identities. We have introduced the 

Hi 
■kG- 

P ^bm 



3H^ 

usual notation pc = p is the baryonic matter and we shall assume 



Pc 



(20) 



The interpretation of equations ( |T71fT9|) is straightforward. Equation ( ITTI) is the Friedman 
equation determining the time evolution of the scale factor a{t). The first term in the right 
hand side of (fT7|) is the contribution from the Born-Infeld field C^^^. Defining the density and 
pressure for the Born-Infeld field, 

_ 1 yM _ 1 XY 

Pbi ~ 7> 7^~v~^^ Pbi — ~- 777^ (21) 
the right hand side of (fTTj) is simply y{pBi + p)- Furthermore, in terms of Pbi and p^j, equation 



(iTSll takes the usual conservation form 

{pBiO?y = -pBiia^y. (22) 

Eq. f|T9|) ("the Friedman equation for the metric q^^") provides the equation of state for pgj 
and Pbj allowing a full solution to the problem. Note that using fl2Tl) the functions X{t),Y{t) 
can be written in terms of pBi{t),PBi{t) and fll9p becomes a (differential) relation between these 
two functions. This equation of state thus have one free parameter represented as an initial 
condition. 

We shall now show that Pbi behaves like dark matter for small times, and like dark energy 
for large times. 

3.3 Asymptotic a ^ and a ^ oo behavior 

Due to the complicated and non-linear character of equations ( ITTlfTQil we shall study them by 
series expansions and numerically. 
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We first study the behavior for large values of a. In this regime, the baryonic matter density 
p ~ does not contribute. (A radiation component would not contribute either.) Neglecting 
the term p/ Pc-, it is direct to see that the functions, 

a{t) = ao e*/^, X{t) = Y{t) = e*/^, (23) 

with C = ^3{1 — a) IHq provides an exact solution to fll7m9p . Thus, de-SitteJ^ space is a 
solution to (I17m9p for large times. The constant C measures the value of the associated vacuum 
density. In order for this solution to approach de-Sitter space with the correct exponent, we 
must impose 

^ ^A. (24) 



3(1 - a)im^ 

Hq and Q\ are determined by observations. This provides a first constraint on the parameters 
/ and a entering in the action. We shall use (l24l) to solve / in terms of a. 

Now, we study the a{t) ^ region. In this regime, an exact solution is not available, but 
one can display a series expansion with the desired properties. The following series 

a{t) = ao t^/^il + 0{t^/^)), X{t) = xl{l + C(t)), Y{t) = Xo(l + 0(t)) (25) 



provide a solution to fll7m9p . The crucial point here is the exponent t^/^ in a{t) meaning that 
by ao- 



C^^p does indeed behave like matter for small times. The amount of dark matter is controlled 



3.4 Numerical interpolation 

Our final goal is to display a solution for a{t) interpolating between a{t) ~ t^^^ for small a{t) 
and a{t) ^ e^* for large a{t). Furthermore, we would like this solution to exhibit the right 
amount of dark matter and dark energy. This will be done by a numerical analysis. 

Equations fll7til9p are of first order and thus we need to give three conditions ai = a(l), 
Xi = X{1) and Yi = 1^(1), plus the values of a and / to integrate them. These are 5 parameters. 
However only two of them are independent. This can be seen as follows. 



^The existence of this exact solution is not at all surprising because we already know that the general 
equations ([5]) and ([6|) accepts solutions of the form i?^^ = Ag^^ when q^^ is proportional to g^^, 
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First of all, for a flat model, we can choose a(l) = 1. Second, in (12^ . we already encounter 
one condition on the parameters to achieve the right evolution. Eq. (IMl) allows to solve / in 
terms of a. One extra condition follows by evaluating Eq. f|T7j) today, 

1 Yl 



1 



hnn 1 



(26) 



from where we can solve X\ in terms of Y\ and /. The remaining parameters are thus a and Y\, 



We have integrated (11711191) numerically varying a and Y\. The resulting curve is compared 
with the evolution predicted by ( |T2lfT3l) . Our conclusions are the following. 

1. First of all, there exists values of a, Yi such that the evolution predicted by f[T^ is almost 
undistinguishable from that following from (I17H19I) . at least for the part of the Universe 
we can observe < t < 1. In Fig. [H the continuous line represents the Born-lnfeld theory 
with a = 0.99 and Yi = 10.59. The dots represents the Friedman evolution dictated by 
([T2|) . The parameters were adjusted for a best fit, and in particular the big-bang occurs 



Born-lnfeld Gravity vs Friedman 




Figure 1: a = 0.99, Yi = 10.59 



aX t = 0.00731.. in both theories. 

2. If a is not close to one, there is no value of Yi to achieve a good fit with the Friedman 
equation. The above picture corresponds to a = 0.99. The value a = 0.9 also gives a 
good fit, but smaller ones do not. a > 1 does not work either. 
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The fact that a ~ 1 to have a good fit is quite pecuhar because the actual value a = 1 
is singular and the de-Sitter solution does not exist (See Sec. 12.21) . In any case, recall 
that a enters in the action as a coupling constant and is not subject to variations. More 
testings on the theory should narrow the actual value of this parameter. 

3. Of course no measurements exist for t > 1, but it is interesting to explore the predictions 
of Born-Infeld theory to larger times. If one chooses the parameters such that the Big- 
Bang occurs at the same value of t in both theories, then for large t the expansion factor 
a{t) grows slightly slower in the Born-Infeld theory. Further details on this issue will be 
reported elsewhere. 

3.5 The evolution of the equation of state 

As we mention in Sec. 13. 2[ the field C^'^p can be characterized by an energy density p^i and 
pressure pgj whose expressions are given in (1211) . The corresponding equation of state is. 



and we observe that the pressure is always negative. Fig. 13.51 shows the evolution < t < 3 
of the quotient Pbi/Pbi- We see clearly the interpolation between p^j = for small times, and 




(27) 
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Figure 2: Evolution of the equation of state 



Pbi = 



Pbi for larger times. 
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4 Spherical symmetry and galactic rotation curves 



In this section we explore the action ([T]) on galactic scales where dark matter also plays an 
important role. (Dark energy is less relevant at this scale.) Stars orbiting galaxies have rotation 
curves not matching the observed luminous matter and this implies the existence of 'dark 
matter'. See [25] for a recent review. If the Born-Infeld theory ([1]) can be regarded as a good 
candidate for this exotic form of matter, it must account for these flat rotation curves. 

Numerical studies of matter interacting only gravitationally suggests the 2-parameter NFW 
[26] density profile for the dark matter halo, 



a 



pNFw{r) = —— — (28) 

The density diverges at the origin but the total mass is finite. The rotation curve associated 

to this profile has a peak and decreases slowly with r. 

Another popular density function is the (pseudo) Isothermic profile, 

h_ 

1 + r^/rp 

having a finite density at the origin. This halo leads to a rotation curve increasing monotonically 
with an asymptotically fiat region. See [27| for an observational comparison of both profiles. 



P'so{r) = T—::nZ2 (29) 



Our proposal is that the Born-Infeld field appearing in ([T]) can account for the dark matter 
present in galaxies. We shall prove that the action ([T]) gives rise to a rotation curve behaving 
like the NFW halo for small r, and as the (pseudo) Isothermal halo for large r. In particular, 
the Born-Infeld theory yields an asymptotically flat rotation curve (at least to first order in the 
coupling constant). To this end, we study in this section solutions with spherical symmetry to 
equations (IMS]). 

4.1 The ansatz and equations of motion 

Our main interest is the 'dark matter' contribution to the gravitational potential induced by 
C^i^p. We shall then neglect here the visible matter, set T'"^'^ = and study the solutions of 



G^,u = -^J%t,aq"'^ gpu (30) 



K^u = ^{9^lu + (yq^u) (31) 
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with spherical symmetry. 

First, note that since f^^jy and q^v are invertible, the right hand side of (l30l) is different from 
zero at all points. This represents a non-compact source. Since for a galactic problem we expect 
the curvature to be very small, the right hand side of fl30l) must be very small. This is indeed 
true because the parameter / is of cosmological scale, and thus is very small compared to 
any galactic scale. 

This also means that we can treat the right hand side of (!30l) and (I3T!) as perturbations. 
The parameter measures the interaction between g^^ and q^y. If — 0, these two fields are 
decoupled and satisfy the order zero equations, 

G,,. = 0, K^, = (32) 

The interactions can be incorporated perturbatively by expansions of the form, 

9^JLy = 5'^°!^ + ^ Q^jll + ji 9^i^l + ■ ■ ■ (33) 

= C + + (34) 

where gj^l and q'j^l satisfy the order zero equations ( l32l) . [There exists a different perturbative 
scheme leading to a different sector of this theory. Instead of treating ^ as very small, one could 
start with the exact solution described in Sec. 12. 2[ and study linearized fluctuations around 
that background. This yields a different set of solutions that will be studied elsewhere.] 

We start our discussion with the order zero fields gj^l and qj^l- The static, spherically 
symmetric solution to the order zero equation Gf^^ = is the Schwarzschild metric 

g°}, dx^dx" = -cM 1 — ]dt^+(l — dr^ + r^dfl^ 35 

However, M represents the total baryonic mass and since we have set T*"^'j^ = we also set 

M = 0. (36) 

Our order zero (7*^^^ metric is thus flat space. The effects of baryonic matter can easily be 
incorporated at the end and will be studied in [20] . 
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The order zero equation for q^j^^ is K^^, = 0. The solution with spherical symmetry is also 
the Schwarzschild metric in the "reciprocal" metric q^iy, 

ql'^l dx'^dx" = ~(3^^ { 1 - ^\ de+{l- ^\ ' ~kV + k\r)dn\ (37) 
V k{r)J \ k{r)J 

Here, k{r) is an arbitrary function of r, Wo is an arbitrary constant with dimensions of length 
(the qfj^iy "Schwarzschild radius"), and /9 is a dimensionless constant. This tensor solves K^^, = 
everywhere except at A; = 0. Some comments on the metric ( |37l) are in order: 

1. We have chosen the radial coordinate r such that gj^l has r'^dQ'^ in the angular part. It is 
then not correct to assume that the metric q'j^l can also be written in terms of r and with 
the same r'^dVL^ in its angular part. This is role of the function k{r) which represents an 
arbitrary radial re-parametrization. This function will be determined by the equations of 
motion. 

2. In the same way, the time scale t is fixed in terms of the metric g^^, and does not need 
to be the same for the metric q^^p. This the role of the dimensionless constant (5 entering 
in ( !37|l . This constant will be important below. 

3. If Wo 7^ 0, the metric ( 1371) solves K^^ = everywhere, except at = 0. As we shall see 
below, the most interesting solution requires wo ^ 0, and explores k{r) all the way to the 
horizon. This means that we are forced to interpret the metric ( l37l) as a black hole in the 
qfj^u space. We shall restrict our discussion to the region, 

k > Wq. (38) 

4. For t^o 7^ it will be convenient to use a dimensionless radial coordinate, 

k = —. (39) 

Wo 

In particular the horizon is now located at, 

k = l, (horizon). (40) 
From now on, all formulas refer to this coordinate. 
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Having chosen the zero order solutions to fl30l) and fl3ip . we now discuss the corrections in- 
duced but the right hand side of these equations. We only discuss here the first order correction 
fo g^u, proportional to p. Since the right hand side of fl30|) is already of order p-, it is enough 
to know q^y to order zero. [Note that qj^l contributes to g'j^^l, qj^l contributes to g'j^l, and so on.] 

Our problem then reduces to replacing g^,^ given by (1371) in (130|) and solve for the metric 
g^j^i, to first order in p. The metric g^i, must be spherically symmetric. We then write, 



-1 



ds' = -c^l + \^r))dt^+(l-'^^j^] rfr^ + rW, (41) 



with 

$ = $(0) + 1 $(!) + ! $(2) + ... (42) 
1 1 

m = m'°^ + — m'^^ + — m<'' H . (43) 

As we have already discussed, in the approximation with no baryonic matter, the zero order 
solution is simply flat space and thus 

= 0, m''^ = 0. (44) 

To first order we obtain the equations, 

dm^- ^ ^^^^ 

(46) 
(47) 



dr 



2(3 dr 


= 0, 


kj dr 


= 0, 


2 

+ fJc wor — 
dr 


= 0, 



dr 

where we have re-defined $^^'(r) in terms of a new function n'^'(r) by 

r— — = u'''(r)H —. (48) 

dr r 

[To first order, the equations only depend on $' and this is why this redefinition does not spoil 
locality.] 
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4.2 Full parametric solution. Two branches 

Equations (I45ti47p are three non-linear equations for the three unknowns m*^^(r), u*^^(r) and 
k{r). A much simpler set of equations can be obtained by changing the independent variable 
from r to k. 

We define the functions u^^\k) , m^^^ (k) and r{k). Also, for any f{r), 

dr dk I dk 

Performing these substitutions, equations fl50tl52p become linear for the unknowns m*^^(/c), u^^\k) 
and r{k), 

/5cW(i-^)^ + 2mW = (51) 

^ + /3c'wor = 0. (52) 

Note in particular that m*^' has decoupled from m*^' and r{k). The general solution can be 
found in closed form, 

r{k) = Ao(-(k-l)\n(l-]-) -l)+Bo(k-l) (53) 

(54) 



27 V kj / ^ V 2 



^{1) ^ + 1^2 ^ + ln(A; - 1) - /lo) (55) 

where Aq, Bq and ho are integration constants. This solution is real for k > 1, that is outside 
the horizon in the reciprocal space 

To explore the properties of the different solutions we first note that the function r{k) 
displayed in (|53|) diverges at two different values of k, 

k = oo, and k = 1. (56) 

Since the function r{k) is a coordinate change and must be globally defined at least in the 
range < r < oo, the derivative dr/dk must be different from zero everywhere. Note that if 
Aq and Bq have the same sign, then the function r{k) has a maximum or minimum, which is 
not allowed. This leaves two simple cases: 
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• Linear branch: Aq < 0,Bq > 0. In this case, r{k) diverges for large k, and becomes 
zero at some finite value k = k^. For k < ko the coordinate r{k) is negative and thus this 
region is not physical. For large k the solution rapidly approaches a linear behavior. The 
physical range of the coordinate k in this case is 

ko<k <oo. (57) 

However, it can be easily seen that m^'^^ /r diverges quadratically for large r. This behavior 
is unacceptable. The same divergency is observed for the potential $(r). From now on 
we shall exclude this case. 

• Logarithmic branch: Aq > 0, i?o < 0. In this case, r{k) diverges at /c = 1 and becomes 
zero at some finite value k = ko. The physical range of the coordinate k in this case is 

l>k>ko (58) 

The most salient and peculiar property of this branch is that infinity is mapped to the 
horizon in the metric There is a strong/weak relationship between both fields. The 
details of this branch are studied in the following paragraphs. 

Fig. [3] shows the behavior of the function r{k) for each branch. 



Linear Branch Logarithmic Branch 




Figure 3: Two branches 
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4.3 The logarithmic branch and asymptotically flat rotation curves 

The most important property of this branch is that the rotation curves are asymptotically flat. 
Let us recall the relation between the Newtonian potential appearing in (HTl) and the rotation 
speed of a (non-relativistic) object at distance r, 



v(r) = \r^^. (59) 
V dr 

(This follows from the geodesic equation.) On the other hand, the derivative of the potential 
$, to first order in p, is given in terms of m'^' and m*^^ in (H8l) . The rotation curve can be 
expressed ClS db parametric function , 



v{k) ^ 



M(^)(fc)+ r = r{k) (60) 



l\ ' ' r{k) 



where u'-^\k) , m'-^^ (k) and r{k) are given in fl53tl55p . 

From these expression is it direct to compute the limit, 

'"Ic = limf^(A;) 

wo{(3^Al-Aw^o) 2 
= AA^PP ' ^''^ 

which is indeed finite. 

However, this is not the whole story. We need to impose boundary conditions at r = 
{k = ko) to ensure that the solution and in particular the rotation curve (!59l) is well-behaved 
there too. This will imply the following constraints and redefinitions of the parameters Aq, Bq 
and Hq. 

1. We first express Bq in terms of ko, the point where r(/co) = 0. This gives the following 
expression for Bq, 

2. Second, must be finite at r = 0. This implies that m'-^\k) must vanish at k = ko and 



this fixes ho to be 



^0 = ^kl + ^kl + ko + Hko - 1) (63) 
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3. Finally, the orbital velocity of an object at r = must be zero. This implies that 
^(1) _|_ 777,(i)/7- evaluated at k = ko must vanish. This is achieved by choosing the constant 
Aq to be 

A, = (64) 

In summary, boundary conditions at r = fix i?o, and Hq in terms of a new parameter k^. 
The full solution is then characterized by three remaining constants. The length scale wq, and 
two dimensionless numbers /3 and kg. 

4.4 A better parametrization and examples 

The solution we have found is still parameterized by several numbers. The functions r{k),v{k) 
depend on /, c, j3, wq, k^. The first two, Z, c enter in the action and cannot be varied. In fact / has 
been already constrained by the cosmological analysis. The other three remaining parameters 
can be chosen to match a desired physical situation. Before plotting examples is it convenient 
to choose a different basis for these three arbitrary parameters. 

First, the asymptotic velocity foo computed in ( 1611) in terms of k^ is 

2 _ 4fco - Akl + k'^-l wl 2 

2{2k,-l)kl ^^^^ 

This parameter is of course a natural observable which can be identified easily for most galaxies. 
We use this equation and express wq in terms of Woo, 



Wo 



^ 4fcg -4kl + k^-l c 
Second, the dimensionless parameter f3, which enter in ( 1371) . can be redefined as 

P=-^. (67) 
ro c 

where Tq is an arbitrary parameter with dimensions of length. 

With these definitions, the functions r{k),v{k) take the convenient form 

r{k) = rofi{k, ko), v{k) = v^f2{k, ko). (68) 

The arbitrary constant tq sets the length scale while v^o set the velocity scale. Since both 
are arbitrary, they can be fixed to any desired values to fit realistic curves. The constant 
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ko controls the shape of the curve and how fast it grows. Since there are three independent 
parameters, there will be a degeneracy when fitting these curves with observational data (this 
will be discussed in |20] ) . The explicit expressions for fi , /2 are not very illuminating, and can 
be derived directly from the solution (15311551) . Of course /2 satisfies /2(1, ko) = 1. 

Fig. (j4]) shows examples of the curve with v^o = 100km/ sec, Tq fixed, and varying /cq- The 
top curve corresponds to ko = 1.5. As ko increases we observe a slower growth of the rotation 
curve. All curves asymptotically reach the value Voo = 100km/ sec. The horizontal axis is 
expressed in terms of r/ro, and choosing tq one can fit any desired length scale. 



lOOn 




" 5 10 15 20 25 30 35 

r/r_0 

Figure 4: Rotation curves for ko = 50, 15, 5, 1.5. 

It is interesting to note that for values of ko smaller than ko ~ 1.5, the curves change shape. 
Fig. ([5]) shows the rotation curve for ko = 1.5, 1.03, 1.003, 1.0005. The top curve corresponds 
to ko = 1.5. As ko becomes smaller, the rotation curves growths more slowly. 

Note that one does not expect the curves to be asymptotically flat to all orders. The 
solutions discussed here are only the first order approximation in the coupling The next 
orders are necessary to extrapolate the result to large values of r, comparable with /. Also, 
the near horizon region for the metric g^i, is singular in Schwarzschild coordinates and thus a 
proper analysis in regular coordinates may also change the behavior near infinity. 

4.5 Final remarks 

We end this section with two extra comments regarding the solutions with spherical symmetry. 
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Figure 5: Rotation curves for ko = 1.5, 1.03, 1.003, 1.0005. 



Orders of magnitude and Solar System: The solutions we have considered contain a length 
scale, wq. This parameter was replaced in (l66l) by the final speed Voo, which is a better 



observable. It is however interesting to estimate the values of wq for a realistic situation. 



We set / ~ lO^kpc (cosmological length), and ^ ' 
Voo ~ lOOkm/sec. Fig. [6] shows wq as a function of ko. 



;10 ^, for a typical situation with 




Figure 6: Wo(kpc) as a function of /cq- 

For ko > 3, Wq is equal to a few kpc. This is a natural galactic scale. With an optimistic 
viewpoint one can thus assign to Wq some physical meaning determined by the length of the 
object observed. In other words, the tensor q^^ is a field whose natural length scale of variation 
is determined by the object. 
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Now, the natural dimensionless parameter which controls the corrections from flat space 
is If we believe that the value of Wq is comparable to the object of study, then for Solar 
System experiments ^ is too small, and the effects of C'^^ should not contribute. 



Central density: The central density associated to C^^^ diverges linearly, as the NFW profile 
(!28|) . This can be seen by solving fl45ti47p . for small values of r, as a series expansion. The 
series, 

k{r) = ko - 7 r + 0(r^) (69) 
mW(r) = -^^iM^ r + C>(r2) (70) 
««(r) = + (71) 

solve fl45tl47p with the boundary condition t'(r) ^ as r — 0. With this solution at hand we 
can compute the behavior of the associated mass density, 

A7rGp{r) = \{r^^y (72) 



r 



P r 

with a linear divergency, as anticipated. 



2ik, - l)^oC^ ^ ^^^^ ^^3^ 



5 Eddington action, the equivalence principle and g^^jj = 

Our proposal for dark matter and dark energy is summarized in the action ([1]). Once the action 
is written one can "roll down" exploring its predictions and consequences by usual methods. 
This is what we have done so far. However, it is also interesting to "climb up" and attempt a 
derivation, or at least a good motivation to include the Born-Infeld term in the gravitational 
action. 

We start this section recalling a well-known effect. Consider a system of spins. If no 
external field is applied (and the temperature is not too small) the macroscopic average is 
{S) = 0. On the contrary, in the presence of an external field, Hext, the symmetry is broken. 
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the spins align and produce a non-zero macroscopic average (S) ^ 7^ 0. It then follows that 
the total magnetic field felt by a charge q is 

HT = He.t + {S)ff^^^. (74) 

The orbit of the charge will obey the Lorentz equation with Ht not Hext- If we did not know 
about spins the contribution (S) ^ would be interpreted as a sort of 'dark' magnetic field. If 
the temperature is below the Curie temperature, the external field could be removed and the 
spins remain in their 'ordered' state with {S)o 7^ 0. 

Let us now describe an analog of this effect in the theory of gravity. Topological manifolds 
are invariant under the full diffeomorphism group. Riemannian manifolds are invariant only 
under the subgroup of isometries of the metric. The state g^^ = represents the unbroken state 
of general relativity [28], and the introduction of a metric breaks the symmetry. The natural 
geometrical analog of the external field Hext is the metric tensor (7^^. (See [291 EQl [31] for other 
discussions on the state g^^i, = 0, and [32] for a recent critical viewpoint.) 

We shall treat the metric as an external field which can be switched on and ofl^. Our first 
goal is to explore fields that can be defined in the absence of a metric. The simplest example 
is given by a connection C'^j^p(x). In fact, Eddington introduced a purely affine theory a long 
time ago [21], 

Io[C] = K Jd'^x sJk,,{C) (75) 

where K^^ is the curvature associated to the connection C^^^(x) (see Eqn. ([2])). This action is 
invariant under spacetime diffeomorphism and yields second order differential equations for the 
field C^'j^p. Eddington action was extensively studied as a purely affine theory of gravity, and 
also as a possible unification of gravity and electromagnetism [2^ [53] . We take here a different 
interpretation and let the field C^^p be an independent degree of freedom. 

We now turn on the external field gfj,iy and study the effects of both g^^i, and C^j.^ on particles. 
The first problem is to determine the action for the coupled system. We do not want to introduce 
ghost or higher derivatives. The action ( 1751) is already free of anomalies. So we start by adding 
the standard Einstein-Hilbert action for g^^^, and consider 

j d\[^R + K^,). (76) 



•^In this picture, the big-bang could be understood as a smooth transition from a manifold without metric 
into a Riemanian manifold. 
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With this action, the fundamental fields g^^i, and C^^^ are decoupled. To make the theory 
more interesting we add interactions. The most attractive theory (although not unique) having 
second order field equations is the Einstein-Born-Infeld action introduced in Eq. ([1]). 

An important point now is to define the geodesic equation for the coupled system. In the 
presence of a metric Qf^u there is a natural affine connection T^^^^ represented by the Christoffel 
symbol, 

r% = ^9^''igau,p + g.p,. - g.p,a). (77) 

The question is, should geodesies be defined with respect to C'^^^, T^^^^, both? In order to comply 
with the equivalence principle we shall postulate that particles only couple to the metric and 
not to the connection C^^^. The geodesic equation then take the usual form 

xf" + T^^i^x'^x^ = 0, (78) 

where T^^^ is the Christoffel symbol ( |77l) . Observe that the metric satisfies the equations ([5]) 
and is coupled to the field C^^p. In this sense, C^^p does contribute to gf^u and indirectly affects 
the motion of particles. This is how the field C can explain fiat rotation curves. 

Now, the analogy with spin systems can be pushed a little bit further. We have seen in 
the cosmological analysis that for large times the system approaches the de-Sitter solution (see 
Sec. 12. 2p . and in particular the metric q^i, becomes proportional to g^^, q^u ^g^u- One 
can interpret this fact as analogous to the alignment of spins along the direction of the applied 
field, {S) XHext- Of course, to support this interpretation one would need to consider generic 
initial conditions. This will be analyzed elsewhere. 

Finally, recall that when the external magnetic field is removed, spins can have a sponta- 
neous non-zero average (S), and this vector generates forces on charged particles. Is there a 
gravitational analogue to this effect? The gravitational force is measured by the connection 
( 1771) . entering in the geodesic equation. The external field is the metric. Now, as the metric 
is removed, the Christoffel connection becomes ^, with the same scaling weight in the numer- 
ator and denominator. For a large class of paths the limit is a finite function. Since the only 
connection available at g^^,^ = is C^^p, it is tempting to conjecture that F^^,^ — > C^^p, as the 
metric is removed. In this way, the geodesic equation has a non-trivial limit when the metric 
vanishes, and particles will feel 'forces'. These are not forces in the usual sense because there 
is no metric. (Although note that a geodesic equation, defined by parallel transport, can be 
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introduced without a metric.) The hmit (7^,^ — > was the key ingredient employed in [22] for a 
different approach to understand dark matter as an effect associated to a topological manifold. 
To make these ideas precise a theory describing the process g^^, — > is necessary. We hope to 
come back to this interpretation in the future. 



6 Conclusions 

Dark matter and dark energy have unique properties and their understanding in one of the most 
crucial problems faced by theoretical physics today. Dark matter does not interact with normal 
matter and this property has motivated us to look for fields which have this property somehow 
"built in" . We have explored gravity coupled to connection C^^^ field with a Born-Infeld action. 

This theory comply with the main background properties normally attributed to dark matter 
and dark energy. First, the evolution of the scale factor in cosmological models has the right 
time dependence interpolating between pressureless matter and a cosmological constant. 

At galactic scales dark energy is less relevant but dark matter still plays an important role. 
By an approximation valid for distances much smaller to the Hubble radius we have solved 
the equations of motion for spherical objects and find the expected rotation curves. These 
curves satisfy the basic asymptotic flatness observed in galaxies providing new support for this 
proposal. 

We have left several topics for the future. The stability of this theory and the study of 
primordial fluctuations are important to determine the CMB anisotropies. This will be reported 
in [19]. On galactic scales a systematic fit with observational curves is necessary. This issue is 
presently under study and will be reported in [20] . 



7 Appendix. Derivation of the equations of motion 



The fields which are varied in the action ([T]) are the metric (yf^j, and the connection F^^^. The 
equations of motion for the metric follow by a straightforward variation of the action. The 
result is 
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(79) 
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This equation can be drastically simplified by using the equation of motion for the connection 
r^j^p. This equation is derived in two steps. First, since the action only depends on the curvature 
K^jj once can compute the variation using the chain rule, 

JL= f ^^(^l^) (80) 

Just like in Eddington [21] theory one finds by direct variation that the combination 

1 r. ZTT^f 1 V 



^^,,p. -^f^-^yg_l2j^J 

satisfies 

Dpi^qn = (82) 

where Dp is the covariant derivative built with the connection Fq. Since F'^,^^ is symmetric, this 
equation imply 

Kp = ^?^"(9a.,p + Qap,. - q.,p,a) (83) 

We thus write F^^^^ in terms of g^^,. The equation (ISTl) now depends only on q^^. Taking the 
determinant at both sides, and inverting one readily derives ([6]). 

The final simplification follows by noticing that the right hand side of (179|) contains ^/qq^'^ ■ 
Thus, using flHT]) . Eq. fl79l) is transformed into ([5]). 

The analysis of these equations is greatly simplified by using the bi-metric formalism 
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